Three-dimensional correlated-fermion phase separation from analysis of the geometric 

mean of the individual susceptibilities 
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A quasi-Gaussian approximation scheme is formulated to study the strongly correlated imbal- 
anced fermions thermodynamics, where the mean-field theory is not applicable. The non-Gaussian 
correlation effects are understood to be captured by the statistical geometric mean of the individual 
susceptibilities. In the three-dimensional unitary fermions ground state, an universal non-linear 
scaling transformation relates the physical chemical potentials with the individual Fermi kinetic 
energies. For the partial polarization phase separation to full polarization, the calculated critical 
polarization ratio is P c = [1 - (1 - 6/5 ]/[! + (1 - C) 6/5 ] = 0.34. The £ = 4/9 defines the ratio of 
the symmetric ground state energy density to that of the ideal fermion gas. 
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Phase separation resulting from the microscopic dy- 
namics involves a wide range of natural phenomena. Cur- 
rently, it also serves as the pivotal topic in the strongly 
interacting fermions quantum many-body theory. 

In the past few years, considerable efforts understand- 
ing the crossover physics from Bardeen-Cooper-Schrieffer 
to Bose-Einstein condensation(BCS-BEC crossover) with 
ultra-cold atomic fermi gases have been made. At the 
Feshbach resonance point, the divergent S- wave scatter- 
ing length with the existence of a zero-energy bound state 
can exhibit the thermodynamic universality [J 

Another exciting issue is on the asymmetric fermions 
phase diagram with unequal populations, which per- 
sists as a fundamental theme for a long time @, i, 0. 
With the Feshbach resonance techniques, tuning inter- 
action strength and controlling the population or mass 
imbalance among the components offer a playground 
to test the non-perturbative many-body theory for the 
asymmetric fermions thermodynamics 0, 0]- The exper- 
iments show that at zero temperature the system can 
undergo a first-order quantum phase transition from a 
fully-paired superfluidity to a partially polarized normal 
gas[8|. Furthermore, when the imbalance population ra- 
tio P = (iVj — Ni)/(Nf +N±) between the two spin com- 
ponents reaches a critical value Pc, the state can transfer 
from the partially polarized phase separation to the fully 
polarized normal one. 

The three-dimensional strongly interacting asymmet- 
ric fermions ground state is one of the prominent un- 
resolved problems^. Meanwhile, the zero-temperature 
ground state energy of the partial polarization phase 
separation is of particular importance for explaining 
experiments ITrjj] . For the equally populated spin states, 
the Monte Carlo methods are believed to be capable of 
doing the simulation calculations. For the asymmetric 
system, the numerical simulations will suffer from the 



serious fermion sign problem Making a definitive 
theoretical conclusion or finding an unbiased analytical 
solution remains as an intriguing but herculean task. As 
far as we are aware, until recently, any analytical attempt 
was not known. 

There can exist a relation between the Pc and dimen- 
sionless constant £ at unitarity with \a\ = oo 
13, 14, HH. The £ defines the ratio of the symmetric 



fermions energy density to that of the ideal noninteract- 
ing ones(T = 0). By exploring the complex non-Gaussian 
correlations, the purpose of this Letter is offering a novel 
tool to calibrate the ground state energy of strongly inter- 
acting asymmetric fermions and fix the relation between 
£ and Pc- The method is based on extending the for- 
malism developed in [ill 17, [ijj]. The calculations are 
performed with ks = ft = 1. 

To address the imbalanced fermions thermodynamics, 
we formulate the universal medium-scaling Hamiltonian 
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In Eq.JT]), a, (3 =| (a), I (b) represent the hyperfine spin 
projection Ising- variables. The Hamiltonian is the same 
as the Bethe-Peierls contact interaction, except that the 
bare coupling constant Uq = 4ira/(2m r ) is substituted 
by a medium regulated functional U*q[ti, T]. The m r — 
m a mb/(m a + nib) is the reduced mass. In the vacuum 
limit n — > 0, the H reduces to the original version pos- 
sessing a global U(l) or Z2 gauge symmetry (assuming 
m a = rat without a loss of generality). 

Due to the medium dependence of U*^ 7 the non-trivial 
87i oc (x r>a [n,T}N a enforce the energy-momentum con- 
servation law[19|. Routinely, the strong correlation ef- 
fects arc further understood as a spontaneously gen- 
erated one-body correlation potential from the density 
functional theory viewpoint [la. llTl] . with which the gen- 
eral thermodynamic consistencies are strictly guaranteed 



2 



within the non-perturbative procedures. 

The medium-scaling formalism allows a natural imple- 
mentation of the strong correlation physics. Initially, we 
need to constitute the U*g functional. The main obser- 
vation is that the environment frustrating characteristic 
can be realized by the twisted contact interaction 

Uo 



U. 



eff 



r - (2) 

With Eq.@, one can define the in-medium inverse scat- 
tering length notation a e g according to U*^ = Aita^/m. 

With the "negative" sign in the denominator of Eq. ([2]) , 
the interaction is alternatively renormalizcd by the strong 
correlation effects or the background fluctuations induced 
in the counterpart Fermi sea. To a great extent, the sus- 
ceptibility Xab is interpreted as an external field mod- 
ulator simulating the mutual non-linear response. Eco- 
nomically, it characterizes the complications of the single 
particle spectrum properties and correlating couplings. 

At unitarity, the scaling thermodynamic quantities de- 
pend solely on the susceptibilities. To quantify the com- 
plicated correlation effects beyond the Gaussian statis- 
tics, the global variable of interest is captured by the 
statistical geometric mean (instead of arithmetic mean) 
of the individual Xa and Xb for the | and J, subsystems 



Xab = y/XaXb- 



(3) 



The susceptibility Xi itself can be separately calcu- 
lated from the Lindhard correlation response function in 
terms of the random phase or s pin "wave" approxima- 
tion, which is also the final goal[18[. Alternatively, it is 
more convenient to calculate the susceptibility Xi with 
the generalized Ward-Identity 



Xi 



(4) 



The fj{z!j) is the Fermi integrals with j = — 5,5,-" 
and the thermodynamical de Broglie wavelength is A = 
-J (2ir)/mT. The collective dynamical variables n* are 
defined by the gap equations Eqs.([7]). The effective fugac- 



ity z: 



1 _ p Pi4 ; 



is analogous to the fugacity z = 1(| 17 1 
At T = 0, Xab is the density of states(DOS) geometric 



Xab = y/N(e a )N(e b ). 



(5) 



Here, N(ei) is the familiar un-perturbated DOS near the 
Fermi surface for the one component fermionsQ. 

From the action Eq.fl]), the grand thermodynamical 
potential 0(T, /i a , /^,) or pressure can be presented as 
the coupled parametric equations with the instantaneous 
quasi-Gaussian approximation method[l6|, [l7| 
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Ma = Ma + 
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n b + fir >a , Mfc = Ma(a ^ b). (7) 



The dynamical gap equations Eqs.([7]) for the corre- 
sponding single particle Green function give the defini- 
tion of the effective chemical potential fi* , respectively. 
The number density n» is expressed in terms of the de- 
fined quasi-particlc Fermi-Dirac distribution function /^^ 



Hi = ^3 /§(*,')> 



(8) 



Sticking to the effective interaction Eq.@ renormal- 
ized by the mixing susceptibility Eq.Q with Eq.([4]), the 
strengths of the correction terms to the Gaussian term 
are self-consistently derived H, TlJ 
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From the underlying grand thermodynamical potential 
Eq.® and with Eqs.([7]), one can derive the remaining 
thermodynamical quantities. For example, the energy 
density e = E/V and entropy density s = S/V read 
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A few remarks are re-emphasized on the grand thermo- 
dynamical potential because it plays an important role. 
The f2(T, /i a , fj, b ) is not the naive polynomial expanded 
according to the bare vacuum interaction strength Uq. 
The collective dynamical variable /x* mixes the low and 
high order contributions. Hidden in the factors C(T,/i*) 
and T>(T,fi*) as well as in Xab, the collective correla- 
tions are combined with the infinite individual dynami- 
cal high order effects. The equations Eq. ©-Eq. ([7]) are 
highly non-linear although they appear a set of coupled 
algebra ones. The fi* dependence of Q can be numerically 
eliminated in favor of the physical chemical potential 

In the following, we will examine the zero-temperature 
asymmetric interacting system thermodynamics, which 
can be derived analytically. At T = 0, the pressure 
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Eq.© and energy density Eq.© expressions read 
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where the "Fermi" kinetic energy e i; particle number den- 
sity n j and rearrangement factor C(cj) are 



2m' 



67T 2 ' 



Cfe) = — . (12) 
The reduced chemical potentials Eq.([7]) and a e Q are 



Ma — £a H «6 + ^(eaj 
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The ki is the corresponding Fermi momentum. 

With unitary a = |oo|, the above expressions can be 
further reduced to an oversimplified compact formalism 



P = T n aHa 
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T n b^b, 
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The salient feature of the main results Eq. (fT3"|) is that 
they appear as the | (a) component contributions "plus" 
those of I (b). The zero-temperature unitary fermions 
thermodynamics obeys the new form of universality, i.e., 
the Dalton partial pressure law of the non-interacting 
ideal gas. Implicitly, the interaction and collective corre- 
lation information is incorporated through the chemical 
potentials Eq. (fT4"|) . 

Let us further discuss the physical chemical potentials. 
They are related with the Fermi kinetic energies of the 
two subsystems through the non-linear transformations 
Eq. (fl"4f . With Eq. (fl"4|) . the energy density and pressure 
can be expressed as the regular functions of the Fermi 
kinetic energies according to Eq. (fT3"|) . 

With the population imbalance polarization ratio P 
(0 < P < 1, assuming n a > n b and /j, a > 



P 



n b 



n b 



(15) 



the numerical solutions of the chemical potentials are 
presented in FigHJ One can see the minority chemical 
potential^) decreases very rapidly while crossing the 
transverse axis. 



FIG. 1: Rescaled chemical potentials versus P at unitarity, 
with "+" indicating the critical position. The above curve is 
for the majority fermions chemical potential while the below 
one is for the minority component. 



For D = 3-dimcnsional, the reciprocal relations of the 
transformations Eq. (|14p are very involved; i.e., it is a 
hard task to express the Fermi kinetic energies in terms 
of the chemical potentials ii a and fi b . In principle, the 
Fermi kinetic energies and particle number densities in 
terms of the physical chemical potentials can be the func- 
tional formalisms such as e a [n a , fJ, b , ( — ) 1 ^ 4 ]. In addition 
to the mathematical complication, tne inverse transfor- 
mation can be singular and nontrivial when fi b — 0. The 
minority fermions chemical potential [i b will change sign 
at this point indicated by Pc- In other words, the minor- 
ity component will quickly "collapse" due to the strongly 
attractive interaction and the system consists of a single 
Fermi surface; below it, there can be two distinct Fermi 
surfaces and the phase separation is favored [2oj|. 

Namely, the vanishing fi b determines the transferring 
criterion position from the phase separation-partial po- 
larization to the full polarization state. The analytical 
relation between Pc and £ can be derived from Eq. (fT4|) 



Pa = 



1_ (1 _fl6/5 
l + (l-£) 6/5 ' 



(16) 



With £ = | , the calculated critical ratio is Pc = 0.34. It 
is in agreement with the Monte Carlo(MC) result Pn » 
0.39 HI- The mean-field theory gives P c « 0.93|23|, 
while the recent experimental value is Pc ~ 0.36§. 

The complete ground state energy versus the polariza- 
tion ratio P is indicated in Fig. [2] and Fig[3l_They are 
in agreement with MC calculations very well[13L l2l|. 

For the unitary fermions ground state, the equation of 
state can be scaled as|l(j, [l2|, [l4[ 



e(n ai n b ) = 
P{Ha,Vb) = 
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with the analytical scaling function 
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g(x) = (1 



1^/6^5/3)3/5^ 



(17) 
(18) 

(19) 
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P c =0.34 



tion and correlation effects beyond the canonical Gaus- 
sian techniques are encoded with the statistical geomet- 
ric mean of the individual susceptibilities. The sim- 
ple medium-renormalized effective action constitutes the 
bridge towards fixing the energy density functional of the 
strong interaction fermions system. 

The non- linear scaling transformation identity Eq. (|14[) 
relates the physical chemical potentials with the Fermi 
kinetic energies. The reassuring association between the 
critical proportional ratio Pc and universal coefficient £ 
agrees well with the Monte Carlo calculations and exper- 
imental measurements. This non-Gaussian correlation 
perspective will stimulate further insights in understand- 
ing the novel quantum phase separation dynamics. 
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FIG. 2: a), Ground state energy versus P with fixed 
density 13;]. b), function h(y) versus y = /it//ia[I3l- The 
dashed curve is for the non-interacting ideal Fermi gas. 




1 \x c =0.49 



(b) 
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FIG. 3: a), e/(|n a e a ) versus £[2l|; b), function g(x) versus 



It is easy to verify that Pq corresponds to the critical 
concentration defined in Refs. flbl 21 1 with xc = 0.49, 
where the scaling g(x) function takes the minimum value 
as indicated in Fig. [3jb). The convex behavior of g(x) co- 
incides with the experimental measurement [Io| . The nu- 
merical solution of the h(y) versus y = fib/na is displayed 
in Figf2]b. These results contribute to understanding the 
realistic trapped system thermodynamics. 

To conclude, the complicated non-Gaussian fluctua- 
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